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Abstract 



We present a new, semi-phenomenological method to extract the conden- 
sate fraction no (T) from two sets of measured responses of ^He with identical 
kinematics at T < Tc and T > Tc. 
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I. INTRODUCTION 



Over the past decades attempts have been made to determine the condensate fraction 
no(T'<) of hquid ^He, which is defined as the asymptotic hmit of the single-atom density 
matrix (SADM) pi(0, s; T<)/p(T<). 

no(T<) = Jimpi(0,s;T<)/p(T<) (1) 

p(T) = pi(0,0;T) is the number density. We shall use the notation 

T< = T <T,;Ty = T > T,; T^ = T<, r> 

In spite of the increasing quality of data, the various experimental results show an un- 
satisfactory spread in hq. The reason may well be that there is no unequivocal, direct link 
between the above definition of no(T<) and experimental information. The quality of an 
extracted condensate fraction thus depends on the accuracy of the approximate isolation of 
the relevant SADM, which contains the condensate fraction. 

Without striving for completeness we briefly review: 

i) computations from first principles; 

ii) results from the comparison of specific theories and data; 

iii) data analysis with minimal theoretical input; 

We then propose and evaluate a novel extraction method which, except for minor as- 
sumptions, is largely model-independent. 

Computations: In principle stochastic methods for a representative sample of a finite 
number of atoms determine the ground state wave functions and thus pi. Those wave 
functions, and thus the SADM, can be computed up to large values of s and (|I|) thus 
directly provide no (see for instance [ In practice the method is accurate only for 

T = K, and from sufficiently large s one extracts n,o(T = 0). Nevertheless, even for T = 0, 
different stochastic methods produce a spread in of, as much as 25%i. 

For finite T a calculation of pi requires a canonical average 
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pi{Q,s-T)/p{T) = Z-\(3) I dri(llk>2 1 cirfc)(ri,rfc|exp[-/3ifA]|n-si„rfe)) (2) 

with Z{(3), the partition function. 

Resuhs for pi(0, s; T 7^ 0) known to us, are for s ^ 7A and for a number of discrete 
around Tc i, but beyond s ~ 2.5A the computed SADM apparently have sizable inaccuracies. 
This reflects on condensate fractions, which have been extracted from s in the range ^ 4-7A 
where a constant value of the SADM is only reached on the averagei. 

Theory: Efforts to obtain the condensate fraction naturally focuss on the extraction of 
the SADM from an observable. A prime example is the cross section for inclusive neutron 
scattering on He, which is proportional to the response (f){q,y). The latter is a function of 
two kinematic variables, e.g. the momentum transfer q and the GRS-West scaling variable, 
y = (M/g)(z/ — g^/2M), with M, the mass of an atom. For fixed q the variable y is a measure 
for the energy loss zv§'i. 

In order to compare with experiment the above response has to be folded into the ex- 
perimental resolution function E. One may then formally write the response and its Fourier 
Transform (FT) in a closed form (see for instance Ref. 

^th,con.^^^ y-T) = J dy' J dy"Fo{y'; T)R{q, y - y'; T)E{q, y' - y"- T) (3a) 
~^th,conv^^^ s-T) = j rfi/exp-*^V*'^'^°™(g, y- T) = Fo{s; T)R{q, s; T)E{q, s; T) (3b) 

The first factor Fo(?/; T) is the asymptotic limit of the response, and depends on the single- 
atom momentum distribution n{p,T). Its Fourier Transform is just the SADM. Final State 
Interactions (FSI) are formally contained in R, which can be expanded in powers of 1/g and 
which depends on higher order density matrices p„. In practice it suffices to retain only the 
dominant FSI oc p2/q- Finally, E is Experimental Resolution (ER) of the measuring devise. 

In the standard approach one computes from dynamics the FSI factor R, which modifies 
the asymptotic response Fq for finite q (see for instance Ref. [ ^,^J^). In addition many- 
body density matrices which weigh all components are required. We focus on Fq, which is 
related to the momentum distribution n{p] T) and which for T< requires modeling 
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n(p,T<) = no(T<)[(27r)='5(p) + /(p;T<)] + A(T<)n"''(p; T<) 

A(T<) = l-no(T<)[l + /(0;T<)] (4) 

In the above parametrization one finds in addition to the macroscopic fraction no(r<) of 
atoms with momentum p = 0, f{p,T^), the fraction of atoms relative to uq with momenta 
p < Pc in the immediate neighborhood of p = oi. Finally, n"°{p, T) above is the normal part 
of the momentum distribution of atoms with p > Pc] A{T) cares for proper normalization. 

Albeit dominant for p 7^ 0, n^"{p,T^) is not directly measurable and is for low T fre- 
quently assumed to be T-independent: ^o°(p; T'<) ^ nQ°{p; T^) = n{p; T^) . For T> it is the 
measurable n{p;T). Equivalently for the SADM 

Pi(0,s;T<)/p(T<) = no(T<) + G'(s;T<) (5a) 
G{s; T<) = no(T<)/(.; T<) + A(T<) [1 + /(O, T^)]pr{s; T^)/p{T^) (5b) 

Above, G is the difference between the SADM and its asymptotic limit uq, which will be 
called the deficit function. By definition lims^ooG'(s, T<) = 0. 

We summarize: Except for T = 0, there is no accurate information on the SADM for large 
s. For finite s an expression for G, e.g. (|5b[), is model-dependent and without knowledge of 



both pi(0, s] T<) and rao(T<) there is no accurate way to reach the deficit function G. 

Data — analysis: One method used in the past relates the integral of the difference 
between responses for T = in the region y ^ to the fraction of atoms having momenta 
in that region. The method is not accurate because of imprecise isolation of FSI and ER 
broadening in (^ (see for instance Ref. [ |[|). 

Of similar nature are cumulant parametrizations of reponses for T<§. Even if the above 
assumption on n"°(p, T<) is made, the data do not permit the extraction of an accurate 
value of no(T<) and pc- Instead one studies fits for varying no(T<)ll2l0. 
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II. CONSISTENCY RELATIONS AND A NOVEL METHOD FOR THE 
EXTRACTION OF THE CONDENSATE FRACTION. 



We now describe a new method for the extraction of the condensate fraction as it appears 
in the S ADM for T< . We try to avoid the described insufficient information on the latter by 
exploiting data on inclusive scattering on two at precisely the same momentum transfer 
q and energy loss z/, or the related scaling variable y. Fairly recent data of that type exist 
for T = 1.6, 2.3 KBIli (MART data). 

We first define from Eq. ( [HED the ratio of the FT of those corresponding data 
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.pi(o,s;r<)/p(r<). 

Eq. (pal) expresses the T-independence of the ER for the MARI data and as a consequence 
the ratio of actual convoluted data are replaced by deconvoluted ones. Those rj^^^ depend 
on q and are complex, because the FT of the data are complex. 

Eq. ( pBD is of a different nature. It assumes that the FSI function R is at most weakly 
T- dependent. Although verified within the accuracy of the data0'0, we cite weak T- 
dependence on the measured pair-distribution function g{r, T)E which, at least in some 
theories, enters the FSI function R @. The price for 'removing' R and the ER is the appear- 
ance of two SADM, instead of the desired pi(0, s; T^) / p{T^). 

Since the precision with which one ultimately extracts the condensate fraction depends 
on the quality of the input, we discuss the latter. We note that t]^^^ in (|^) uses not actual 
data, but instead their FT. With substantial noise, in particular in the larger y tails of the 
data for (j){q,y), one cannot avoid smoothing, or even cutting out those low-intensity data, 
which affect the accuracy of rj'^'^^. 

We first apply the method to the above-mentioned MARI data at T< = 1.6K;T> = 
2.3 Ki from which we chose a representative sample of highest quality for which q = 
17, 21, 23, 29A ^ . From those we extract the ratio lm[ri^^P{q, s)]/'Re[ri'^^P{q, s)], which is well 



determined in the range s < 2.5A. For growing s, the imaginary part increases from but 
is, even at medium s ~ 2.5A, only 2% of the real part, and can thus for all purposes be 
neglected. This confirms that at least the ratio of the involved FSI functions R is virtually 
T-independent. 

Next we investigate the g-dependence of Re[77'^^^] and show in Fig. 1 data for some 
individual q. Before judging the quality we mention that, contrary to predictions, some 
data for T> and varying q, and more for T<, are not uniformly smooth in q. This may cause 
apparent g-dependence in r] and justifies the use of an average for s < 2.5A 
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Fig. 1. The ratio r/^^'P(s;r< = 1.6K,T> = 2.3 K,) Eq. (||), for q = 
17,21,23,29A" and its average @ over q (drawn line). The long dashes 
give ry^°™P, Eq. (|). 
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As a first step we test (|6BD with the only available computed results on SADM for T^i. 
Without any dynamical cause, the outcome for T=1.54 K closest to T=1.6 K of the MARI 
data, is not very close to values from interpolation between neighboring T. We suspect 
the mentioned inaccuracies which are also apparent in non-vanishing condensate fractions 
for T> and the non-smooth dependence of no(T'<) as function of T outside their estimated 
theoretical error bars. Nevertheless we use interpolated computed results and define for 
T< = 1.6K;r> = 2.3 K 

-pi(0,s;2.3)/p(2.3)i-'"f 



r/'=°'"P(s;1.6,2.3) 



(8) 



Pi(0,s;1.6)Ml.6 

The above as function of s < 2.5A is also entered in Fig. 1 and is seen to represent the data 
rather poorly. We shall thus exclusively use jj^^P'"-'" and return to Eq. ([5a|), which by means 
of and (|BD becomes 

"«(T<)^MM|Z)_G„;T,) (9) 

We note that, contrary to (|I|), Eq. (|D is not applicable in the asymptotic region, because 
in order to produce a constant, numerator and denominator in the ratio in (y) have to tend 
to zero in exactly the same fashion. This is an impossibly stringent demand. 

Although Eq. provides in principle a value of the condensate fraction for 

arbitrary, finite s, it is clearly of greater interest to consider a range of s. In fact, the 
right hand side of @ is a function of s, whereas the left hand side is a constant, implying a 
consistency test, to be passed as a condition for a succesful extraction of no(T<). The size 
of the s-range depends on the available information on the various functions in (P). 

The most delicate source of information is the deficit function G. After the warnings 
above, one clearly does not want to use (^) and has thus to rely on the model- independent, 
not too accurate results of Ref. [ ^, which for 4 < s(in A) < 7, provide nQ{T) with uncertain- 
ties of the order of 15%. The crucial point is that those uncertainties decrease relatively to 
G'(s;T<) which increases for decreasing s. We thus conclude that, up to medium s < 2.5A, 
one can trust and use the computed G. The above happens to also be the range for which 
the available data determine r] sufficiently well. 



We thus perform the consistency test, imphed by (|D and it comes as a surprise that the 
right hand side in the considered s-range is only weakly s-dependent, leading to 

no(1.6) = 0.0625 ±0.0017 

no(0) = 0.090 ±0.030 (10) 

The latter is the result for a standard extrapolation to T = 0. The small error limits are due 
to averaging results for different s; we could not estimate the same for uncertainties in G. 
We note that the errors on the three functions in (P) are uncorrelated, which underscores 
the feducity of the extracted value of no(T<). 

The above numbers use G{s;T = 1.6K), interpolated between values for T= 1.54,1.82K 
as reported by Ceperley. We mentioned that the former set, does not interpolate smoothly 
between values for T=1.18, 1.82 K. Smoothing leads to a lower condensate fraction. Obvi- 
ously the accuracy with which one can extract the condensate fraction depends on the same 
for the input. 

Above we also listed older Argonne data for fixed q = 23. lA ^ and some 10 temperatures 
around T^. Unfortunately we do not have available the actual data and the ER functions, 
and have to analyze instead parametrizations of the above, as given by the authors. Those 
are given by the authors in the form of double Gaussians for data deconvoluted from ER, 
and in addition also deconvoluted from FSI. The information no doubt reduces the required 
accuracy. Yet, in view of scarcity of information we performed the above analysis for the 
Argonne data. In order to compare with the MARl data and we fix T> at 2.3 K. 

One first observes that the computed FT of the above parametrized data are essentially 
real. Next, ratios of those FT for and T< and the above T> are identical within 1%, in 
support of (^. Nevertheless, only a limited part of the information can be used for an 
analysis of the above type. Thus for > T > 1.8 K, uqIT^) is the difference of nearly 
equal terms in (^. A reliable determination requires a precision on pi(0, s; T>)/p(T>) and 
G'(s;T<), which is beyond the Ceperley results. A different difficulty occurs for T < 1 K, 
where one has to make an uncertain extrapolation, using the lowest T results of Ceperley. 

8 



The following outcome 



no(T= 1.0 K) 



0.063 ±0.006 



no{T = 1.5 K) 



0.060 ±0.008 



no{T = 1.8 K) 



0.050 ±0.005 



(11) 



carries some 20-25% uncertainties, which appreciably exceed the one in ([To|) based on the 
MART data. One notes that r;,o(T ^ 1.6 K) from the MARI and Argonne data approximately 
agree. However, in view of the fact that we had to use parametrizations of the Argonne data, 
we do not attach too much significance to the correct trend of no{T) as function of T and 
the somewhat low, but not unreasonable, average {no{T = 0))^^ =0.079. 

In conclusion, we have suggested and worked out a new method to extract the condensate 
fraction in ^He from the Fourier transforms of data sets on structure functions at the same 
momentum and energy transform, but for two T below and above Tc. In addition the 
method requires some previously computed dynamical information. The expression for 
?T,o(T<) appears in principle as a function of s, but is in practice a well-determined and thus 
meaningful constant, the condensate fraction. 
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